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We argue that the longitudinal lepton polarization and the forward-backward asymmetries in 
exclusive B — » (K,K*)£ + £~ decays are largely unaffected by the theoretical uncertainties in the 
long-distance contributions both induced by the long-range charming penguins in the £ + £~ channel 
and those in the meson channels encoded in the meson transition form factors, and are mainly 
determined by the short-distance Wilson coefficients at the low-energy scale /i ~ mt. Thus, the 
above mentioned asymmetries provide a powerful probe of the Standard Model and its extensions. 

The investigation of rare semilcptonic decays of the B meson induced by the flavour-changing neutral current 
transition b — > s provides an important test of the Standard Model (SM) and its possible extensions. Rare decays 
are forbidden at tree level and occur at the lowest order only through one-loop diagrams. This fact opens the 
y—{ , possibility to probe at comparatively low energies the structure of the electroweak theory at large mass scales, due 
to the contributions of virtual particles in the loops. Moreover, rare b — > s transitions are expected to be sensitive 
& \ to possible new interactions. These interactions govern the structure of the operators and the corresponding Wilson 
coefficients, which appear in the AB = 1 QCD-improved effective Hamiltonian describing the b — > s transitions at 
low energies. 

However, any reliable extraction of the perturbative (short-distance) effects encoded in the Wilson coefficients of 
the effective Hamiltonian Jl|-Q| requires an accurate separation of the nonperturbative (long-distance) contributions, 
which therefore should be known with high accuracy. The theoretical investigation of these contributions encounters 
the problem of describing the hadron structure, which provides the main uncertainty in the predictions of exclusive 
<**T} | rare decays. 

In rare exclusive semileptonic decays induced by the quark transition b — > s£ + £~ one faces the two types of the 
long-distance (LD) contributions. First, these are the LD effects in the £ + £~ channel generated by the long-range 
charming penguins, i.e. intermediate cc states described by the cc resonances (ip,tp', ...) and the charmed hadronic 
continuum states. And second, these are the LD effects in the meson transition amplitude of the Effective Hamiltonian 
encoded in the meson transition form factors. 

In this letter we study the influence of the above mentioned two types of the LD effects on observables in exclusive 
*-£h j B — > K(*H + £~ decays. We show that the asymmetries of the lepton distributions remain largely unaffected by the 
uncertainties in the LD effects but turn out sensitive to the values of the Wilson coefficients thus providing a possibility 
to reliably separate and study the short-distance (SD) effects. 

1. Neglecting the strange quark mass the effective Hamiltonian describing the b — > s£ + £~ transition has the following 
*^ ' structure I2B 
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where At = V t * Vtb and sl = s(l + 75). The quantities Ci(fx) are the SD Wilson coefficients which are obtained by 
integrating out the heavy particles. The values of the Wilson coefficients at the scale \x ~ 5 GeV are Q: Ci(/x) = 0.241, 
C 2 (vt) = -1.1, C 7l (fi) = 0.312, C 9V (/i) = -4.21 and CW(/i) = 4.64. 

An important contribution to the b — ► s£ + £~~ transition is given by the intermediate cc states in the £ + £~ channel, 
generated by the 4-quark operators in the effective Hamiltonian describing the b — > s transition, and proceeding 
through the chain of transitions b — > sec — > S7* — * s£ + £~ . The cc pairs in this channel generate both the SD and LD 
contributions. 

The b — > sj* transition vertex is transverse with respect to the photon momentum because of the gauge invariance. 
In the limit m s = within a theory with a V — A charged current it has the general structure of the form 

(5(^)1^(0)16^)) = | e |^A t (Ci + C 2 /3)s L (fc 2 ) [( 7 ^ 2 - qrf)Gitf) ~ 2zm b a^G 2 (q 2 )} b(h), (2) 

where G\ and G 2 are the independent form factors. It is important to point out that the form factor Gi(<7 2 ) has 
no pole at q 2 = 0. This property is related to the fact that the photon mass remains unrenormalized at all orders 
in perturbation theory. Also notice that the pole in G\{q 2 ) would have yielded a nonzero amplitude of the radiative 
B — > if 7 decay. 
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The amplitude of the b — > s£ + £ transition through the cc penguin then has the form 

A(b -» sec -» s 7 * - s^+r) = -^|A t (Ci + G 2 /3)^e 2 s L [l^Gf (« 2 ) - 2im h a^q v Gf (q 2 )] b ■ l lfi l, (3) 

where G\ c 2 are the cc contributions to Gi,2- 

The Lorentz structure of this amplitude is similar to the structure of the effective Hamiltonian (Q), so the cc 
contributions to the b — > s£ + £~ transition can be translated into additions to the Wilson coefficients as follows 



16n 2 

T 

C 7l -> C#V) = G 77 + AG 7 c 7 ( g 2 ), AGf^ 2 ) = -^-(Ci + G 2 /3)Gf (<? 2 ). 



G 9y ^G 9 e y V) = G 9 y + AG 9 c y (<z 2 ), AC 9 t(<z 2 )^^(C 1 + C 2 /3)Gf(g 2 ) (4) 



In the g 2 -region of the vector cc resonances (V) the form factors G\ c 2 have singularities corresponding to these 
resonances. The general expression for the b — > sV amplitude can be written as (cf. 

A(b - aV) = -^|A t (Ci + G 2 /3)s L [ 7p M 2 5l (Af 2 ) - 2^m b CT Mt ,<f g 2 (Af 2 )] be^(q), (5) 

where e^(g) is the V polarization vector. Then the form factors GJ c 2 (<z 2 ) in the region q 2 ~ My are expressed through 
the <?i,2 as follows 

G ? S (« 2 ) = 9 i M ^l , f /vM v + [regular terms at g 2 = M 2 ]. (6) 

My — q^ — ll v My 

The leptonic decay constant fy is defined by the relation (0|c7 M c|V) = eY^(q)Myfy. 

If we neglect the soft nonfactorizable gluon exchanges, the amplitude of the b — > sl + l~ transition through the 
charming penguin is connected with the charm contribution to the vacuum polarization, , defined as 

KM = 7 / dxe-^(0\T(c(QhAO),c(xh,c(x))\0). (7) 



i 

As a consequence of the vector current conservation, IK° has the transverse structure 

n£ = (9^q 2 - q^uW%q 2 ). (8) 

The quantity II cc (g 2 ) contains both the LD and SD contributions and does not have a pole at q 2 = because of the 
gauge invariance. In this factorization approximation one has 

Gf = n cE (g 2 ), Gf = 0. (9) 

Similarly, assuming factorization for the b — > sV amplitude, one finds 

gi(M 2 ) = f v /M v , 52 (M 2 ) = 0, (10) 

and relation (^|) takes a simple form 

°f = T72 ifc T , f + [ re S ular terms at « 2 = M vL G ftf) = °> (11) 

My — q z — ll V M V 

where the regular terms come from the contribution of other resonances and the cc hadronic continuum states. A 
standard procedure |^|,^| relies in assuming the cc hadronic continuum to be described by the c— quark loop in the 
spirit of the quark-hadron duality. Then, summing over all relevant resonances and re-expressing the decay constants 
through the corresponding leptonic decay rates with the relation 

r(v^/+r) = 7ra 2 16 f v 



' 27 M v 

we come to the following representation for the shift of the Wilson coefficient C 9 y in the factorization approximation 
(cf- |I) 
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AC^(<z 2 ) = [3Ci(ro t ) + C 2 (m b )} ■ 

The fudge factor k is introduced in Eq. ( |l2| ) to account for inadequacies of the naive factorization framework (see 
H for more details). Phenomenological analyses || suggest that in order to reproduce correctly the branching ratio 
BR(B -> J/ipX -> = BR(B -> J/ipX) • BR(J/^ -> it should satisfy an approximate relation 

k [3Ci(mf,) +C2(to&)] ~ — 1. The SD contributions are contained in the function h(m c /mb,q 2 /m 2 ), which describes 
the one-loop matrix element of the four-quark operators (see, e.g., [Q for its explicit expression). 

It is known that the charm contribution described by (|l^) yields a strong interference between the LD and SD 
contributions to the decay amplitude in a broad range of q 2 . It has been argued however that a simple parametrization 
of the hadronic cc continuum by the quark loop in the spirit of the quark-hadron duality considerably overestimates 
the net effect of the charm at small q 2 . One possibility to describe the cc continuum contribution at small q 2 in 
a more realistic way is to take into account the regular terms in (^) by representing the LD contributions to the 
Cgv as a sum of the Breit-Wigner resonances but assuming a g 2 -dependent f%p{q 2 ) @. Another possibility has been 
implemented in Ref. |^| where the cc contribution to Cgv has been described as a sum of the Breit-Wigner resonances 
and the continuum contribution which has been expressed through the observable cross-section of the process £ + £~ — > 
charmed hadrons. 

Since at present there are no firm arguments to determine with a sufficient accuracy the LD cc effects we consider 
the difference provided by the three models of the AGgy, namely, the one given by Eq. jl^), and those of Refs. 
as a typical present-day theoretical uncertainty in Cgv- The corresponding Cg(/ (q 2 ) are plotted in Fig. 1. We point 
out once more that the pole at q 2 = in Cgy which has been discussed in || is ruled out by the gauge invariance. 

Let us now briefly discuss the nonfactorizable effects. First, notice that factorization yields vanishing of the cc 
contribution to the amplitude of the radiative b — > decay at q 2 = 0: in this case G2 = and the structure 
proportional to G\ does not contribute at q 2 = 0. Hence, the LD cc contribution to the radiative b — > transition is 
a purely nonfactorizable effect. It should be also taken into account that a nontrivial cc contributions to the radiative 
b — > S7 decay AG^(O) ^ implies a nonzero AC^(q 2 ) in the semileptonic b — > s£ + £~ decay. 

A stringent test of the factorization or its breakdown is provided for instance by the polarization of ip in the 
inclusive and exclusive b — > sip transitions: the ratio of the longitudinally-to-transversely polarized ip does not depend 
on the Wilson coefficients but is rather a function of g\(M 2 ) and g2(M 2 ). The analysis of Ref. || favours the value 
\92(M 2 ) I gi(M 2 )\ ~ O.l-f-0.2, although the factorization prediction gi{M^) = is also not ruled out. For estimating at 
all kinematically accessible q 2 the nonfactorizable effects which lead to a nonzero ACy° it seems reasonable to assume 
that the ratio G2/G1, which governs the size of these effects, is limited by the relation |G2(<7 2 )/Gi(g 2 )| < 0.2 in the 
whole kinematically accessible region of q 2 . In this case, the size of the nonfactorizable effects in the region of the 
resonances corresponds to the analysis of j^j, and at q 2 — this prescription yields a nonfactorizable contribution to 
the b — > S7 decay within 5% in accordance with the estimates of Ref. p0| . In the following analysis of the sensitivity 
of the asymmetries in exclusive B — > K^*H + £~ decay to the uncertainties in the LD effects we also allow for a 
nonfactorizable contribution |AG^(g 2 )| < 0.2|AG§f (g 2 )|. 

2. The amplitude of the B -> K { *H + £- decay is given by the relevant mesonic matrix element of the effective 
Hamiltonian. The LD effects connected with the meson formation in the initial and final qq channels are encoded in 
the meson transition form factors of the bilinear quark currents from the effective Hamiltonian (|l|) . Various theoretical 
frameworks have been applied to the description of meson transition form factors: among them are constituent quark 
models ||, QCD sum rules lattice QCD @@' analytical constraints ||]. 

Lattice QCD simulations, because of its most direct connection with QCD, are expected to provide the most 
reliable results. However, at present the lattice calculations do not provide the form factors in the whole accessible 
kinematical decay region as the daughter light quark produced in b decay cannot move fast enough on the lattice and 
one is therefore limited to the region of not very large recoils. For obtaining form factors in the whole kinematical 
decay region one can use extrapolation procedures based on some parametrizations of the form factors. For instance, 
in p8| a simple lattice-constrained parametrization based on the constituent quark picture Jl2| and pole dominance 
is developed. 

QCD sum rules give complementary information on the form factors as they can calculate the latter at not very 
large momentum transfers. However in practice various versions of the QCD sum rules give remarkably different 
predictions, being strongly dependent on the technical subtleties of the particular version p4|-^6[. 

Constituent quark models (QM) have proved to be a fruitful phenomenological method for the description of heavy 
meson transitions. A long-standing shortcoming of the quark model predictions for the form factors has been a strong 
dependence of the results on the QM parameters jl^j. However, as we have found recently ^0|, a combination of 
the results of the lattice simulations at large q 2 and the dispersion approach based on the constituent quark picture 
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pl[ allows one to considerably increase the accuracy of the predictions. Namely, constraining the parameters of 
the quark model by the requirement that the form factors of the QM at small recoils reproduce the results of the 
lattice simulations considerably decreases the uncertainty in the relevant QM parameters. Once these parameters 
are fixed, the spectral representations of our dispersion QM allows a direct calculation of the form factors in the 
whole kinematically accessible region of q 2 . The form factors of the dispersion QM develop the correct heavy-quark 
expansion at leading and next-to-leading 1/toq orders in accordance with QCD for the transitions between heavy 
quarks; for the heavy-to-light transition the form factors of the dispersion QM satisfy the relations between the form 
factors of vector, axial-vector, and tensor currents valid at small recoil. In addition they satisfy the dispersive bounds 
and thus, as well as the form factor parametrizations of Ref. ]l8|], satisfy all known rigorous theoretical constraints. 
However some uncertainties are still present which are connected with the errors in the results of the lattice simulations 
at small recoil and the approximate character of the extrapolating formulas used in |18| as well as the constituent 
quark picture used for the description of the form factors in Ref. [g(J. We consider that the difference between the 
form factors of Refs. jl^,^0| provide a typical theoretical uncertainty of our understanding of the LD effects encoded 
in the transition form factors. 

3. We evaluate the forward-backward (Afb) and the longitudinal lepton polarization (Pl) asymmetries for various 
parametrizations of the LD cc contributions to the Wislon coefficients Cg(/ and CV 7 and various sets of the form factors 
with formulas given in p0[] . Fig. 2 shows these quantities evaluated with the GI-OGE set of the form factors from 
|pOf and AC^ = and various prescriptions of Cg(/ (q 2 ). Fig. 3 demonstrates a sensitivity of the asymmetries to the 
nonfactorizable effects parametrized as additions to the Wilson coefficient Cj 1 . One can see that the uncertainties due 
to the LD cc effects in both CV 7 and Cgy provide a minor influence on the observable quantities in a broad kinematical 
region beyond the resonances. Notice that within the SM Pl(B — > if/i + /i~) ~ 2CqvC\qa/ [Cg V + C 2 0A ) ~ — 1 
independently of the prescriptions chosen for the LD effects in the Wilson coefficients. 

Fig. 4 presents the same asymmetries evaluated with the Cgy given by the Eq. (|l2|) and Cy 7 = and the two sets 
of the form factors: namely, the lattice-constrained parametrization of Ref. and the GI-OGE set from Ref. [ p0[ . 
One can observe that the asymmetries turn out to be more sensitive to the variations of the transition form factors 
with respect to the uncertainties in the LD effects induced by the cc penguins. Nevertheless, there is still a broad 
range of momentum transfers in which the interference between the LD and SD contributions is negligible. Thus, a 
study of the asymmetries under consideration in these regions provides a potential possibility to measure the Wilson 
coefficients. 

In fact, a sensitivity of the asymmetries to the SD contributions encoded in the Wilson coefficients might be 
observed. The CLEO data on the radiative exclusive and inclusive b — > s-f decays allow i?7 7 = Cj~,{Mw) /C^ 1 (Mw) 
in the ranges (see, e.g. ^2j) 0.4 < i?7 7 < 1.2, and —4.2 < i?7 7 < —2.4. As discussed in Refs. J2^j2^] scanning over 
the MSSM parameter space correspondent to the allowed C7 7 regions provides only minor changes in Cgv{Mw) and 
Cioa(M\y)- Thus for illustrating possible new physics effects in lepton asymmetries we fix Cgv(M\y) and Cioa(Mw) 
to their SM values, and vary C7 7 in the allowed regions. The corresponding Afb and Pl are shown in Fig. 5. One 
finds the shape of Aps to be sensitive to the C7 7 value (or at least to the C7 7 sign), and this effect far overwhelms 
the uncertainties in the LD contributions. 

A weaker but still visible sensitivity of Pl{B — > K*£ + £r) to the SD coefficients can be seen (Fig. 5. b). 

The longitudinal lepton polarization asymmetry Pl{B — > K[i + pT) is equal to Pl ~ 2CgvCioA/{Cg V + C 2 0A ) if 
C*7 7 lies in the allowed regions at all kinematically accessible q 2 , except for the end-points and regions near ip and ip' . 
Therefore, the measurement of Pl(B> —> K[i + n~) can provide us direct information on the ratio Cgy /Cioa- 

In conclusion, we have analyzed the sensitivity of the lepton asymmetries in rare exclusive b — > s£ + £" transitions 
to the LD effects. We have found the most important uncertainty to come from the errors in the meson transition 
form factors whereas the uncertainties in the LD effects induced by the long-range cc penguins are of relatively minor 
importance. So for obtaining better predictions more accurate calculations of the transition form factors in the whole 
accessible q 2 region are necessary. Nevertheless, there are broad q 2 -intervals where the interference between the LD 
and SD effects is negligible and the asymmetries in these regions contain information on the SD effects. 

In particular, the shape of the forward-backward asymmetry in B — > K*/i + /i~ turns out to be specific for the 
particular values of the SD Wilson coefficients: namely, within the SM, Afb is positive at small q 2 , has a zero at 
q 2 ~ O.I5Af| and then becomes negative at larger q 2 . At the same time, the present experimental restrictions allow 
a region of the C7 7 with a sign opposite to the SM value which yields an essentially different shape of the Afb (cf. 
also p4|). The longitudinal lepton polarization asymmetry Pl{B — > K*fi + n~) also contains important information 
on the Wilson coefficients at the scale fi ~ mi,, and Pl{B — » Kfi + directly measures Cgy/CioA at \i ~ m^. 

Thus, the experimental study of the forward-backward asymmetry and the longitudinal lepton polarization asym- 
metry potentially provides an effective test of the Standard Model and its possible extensions. 

We would like to thank M. Misiak for helpful comment on the MSSM parameter space. 
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FIG._1. The effective Wilson_coefficients Cg(- J (g 2 ) for various parametrizations of the cc contributions: solid - 

The bold solid line corresponds to Cgv without the LD contributions, namely, 
the horizontal line corresponds to ~Cgv(/J. = 5 Gev) = 4.21. 



dotted 



dashed 



Eq. ((TH) 

CW(aO + (3Ci + C2)h(m c /m b , q 2 /mf) 
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FIG. 2. A sensitivity of the A FB in B -> K*£ + £~ (a), P L in B -» K* 



in Cg{/. The GI-OGE Set of the form factors 



and ACy^ = are used 



£ + £~ (b), and Pl in _B — > K£ + £~ (c) to the uncertainty 
The curves correspond to various parametrizations 



(q ). Notations see Fig. 1. 
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FIG. 3. A sensitivity of the A FB in B —* K*l+l~ (a), Pl in B — > K*£ + e~ (b), and P L in 5 ^ Kt^~~ (c) to the 
nonfactorizable effects parametrized as additions to CV 7 : The Cgy from Q and the GI-OGE set of the form factors 
are used. The curves correspond to AC% = 0.2 AC^ (dashed) AC^ = (solid) AC% = -0.2Aq^ (dotted). Bold line 
corresponds to CW without the LD contributions and AC^ = 0. 
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FIG. 4. A sensitivity of the Aps (a) and Pj, (b) in B — > K*£ + £~ to the meson transition form factors: solid - GI-OGE Set 
from Ref. | p0[ ; dotted - lattice-constrained parametrizations from Ref. p^ . The Cq(/ from M and AC?^ are used. Bold lines 
correspond to the Cgv without LD contributions. 
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FIG. 5. A sensitivity of the Afb (a) and Pl (b) in B — > K*£ + £~ to the Wilson coefficient C7 7 : the asymmetries are 
evaluated for the GI-OGE Set of the form factors feci] and different values of Rj~, = CV 7 (Afiy) /Cfj 1 (Mw) from the allowed 
region J24I : upper dotted line - i?7 7 = 1.2, lower dotted line - J?7 7 = 0.4, upper dashed line - J?7 7 = —2.4, lower dashed line - 
i?7 7 = -4.2, solid line - R 7 ~, = 1 (SM). 
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